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We investigate the dynamical chiral symmetry breaking and its restoration at finite density and 
temperature within the two-flavor Nambu-Jona-Lasinio model, and mainly focus on the critical 
behaviors near the critical end point (CEP) and tricritical point (TCP) of quantum chromodynamics. 

The multi-solution region of the Nambu and Wigner ones is determined in the phase diagram for 
the massive and massless current quark, respectively. We use the various susceptibilities to locate 
the CEP/TCP and then extract the critical exponents near them. Our calculations reveal that the 
various susceptibilities share the same critical behaviors for the physical current quark mass, while 
they show different features in the chiral limit. 


PACS Numbers: 12.38.Mh, 12.39.-x, 25.75.Nq 

I. INTRODUCTION 

Quantum Chromodynamics (QCD) is often viewed as 
the basic theory of strong interactions, whose degrees 
of freedom include quarks and gluons. As a non-Abelian 
gauge field theory, QCD exhibits two important features: 
asymptotic freedom and color confinement. It works 
well in the large-momentum-transfer processes through 
the perturbative techniques due to the feature of asymp¬ 
totic freedom. However, in the small-momentum-transfer 
processes, the coupling constant becomes so strong that 
problems have to be treated by many non-perturbative 
methods, such as lattice QCD [1], Nambu-Jona-Lasinio 
(NJL) model [2, 3], Dyson-Schwinger equations [4] and 
so on. 

It is believed that the strongly interacting matter un¬ 
dergoes some phase transition and the chiral symme¬ 
try is partially restored at high temperature (T) and/or 
quark chemical potential (/i), from the hadronic matter 
to the quark-gluon plasma (QGP). This phase transition 
is expected to be produced in on-going heavy-ion colli¬ 
sion experiments, such as the BNL Relativistic Heavy- 
Ion Collider (RHIC) and the CERN Large Hadron Col¬ 
lider (LHC) [5-7]. In addition, high density QCD matter 
at low temperature is anticipated in the compact stars 
physics [8-10]. 

Lattice QCD works well for low {i and finite T. How¬ 
ever, the lattice simulations have trouble in studying the 
QCD phase transitions with finite // due to the severe 
fermion sign problem. Therefore simpler and mathemat¬ 
ically tractable models which respect the essential sym¬ 
metries of the QCD are necessary. NJL model captures 
fundamental features of QCD itself, such as the dynam¬ 
ical chiral symmetry breaking in the vacuum, and can 
therefore yield profound insight into the critical behav¬ 
iors associated with chiral symmetry. 
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Many effective theories [11-13] of QCD predict the ex¬ 
istence of the tricritical point (TCP) and the critical end 
point (CEP) in the QCD T — fi phase diagram in mass¬ 
less and massive quark cases, respectively. TCP is a point 
where three coexisting phases become identical simulta¬ 
neously. The chiral phase transition is of first-order at 
high /i and low T, then turns to be second-order after 
the TCP. The second-order transition line is expected 
to be in the universality class of three dimensional 0(4) 
symmetric spin model, also called 0(4) line. With the 
current quark masses increasing from zero to the physi¬ 
cal masses, TCP emerges to CEP at which the first-order 
phase transition line ends from higher /i toward lower ji. 
Meanwhile, the second-order phase transition turns into 
a smooth crossover at lower fi and higher T in the phase 
diagram. Different from the 0(4) line, the CEP is ex¬ 
pected to be in the Z 2 universality class [14], this is just 
one of several possible scenarios for QCD. 

Both at the TCP and CEP, second-order phase tran¬ 
sitions are expected to occur. Therefore, the values of 
critical exponents will be an interesting and important 
question. The critical exponents in mean field approxi¬ 
mation have been predicted, by expanding the Cinzburg- 
Landau thermodynamic potential to the order parame¬ 
ter [15]. Therefore it is important to verify these uni¬ 
versality arguments from the study of specific effective 
theory of QCD, for this purpose we investigate the crit¬ 
ical behaviors of the QCD by using the NJL model for 
two quark flavors. 

The rest of this paper is organized as follows: In Sec. 
H a mean field description of the Nambu-Jona-Lasinio 
model with two quark flavors as an effective realization 
of the low-energy sector of QCD is presented. In Sec. 
HI we display the multi-solution region of the Nambu 
and Wigner ones and analyze the behaviors of the ef¬ 
fective quark mass and the quark number density in the 
T — {1 plane for massive and massless current quark, re¬ 
spectively. In Sec. IV we analyze the behaviors of some 
susceptibilities in the T — fi plane and determine the lo¬ 
cation of the CEP (or TCP). Then critical exponents on 
the CEP are calculated and compared with those of the 
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TCP. Finally, in Sec. V we will summarize our results 
and give the conclusions. 


II. THE NAMBU-JONA-LASINIO MODEL 

The Lagrangian of the NJL model [2, 3, 16] is written 
as: 


the real time [17]. So the Green function of a free fermion 
at temperature T and chemical potential (i is given as 

G{p; T, fi) ={fi + M){ ^ / . + 27ri6{p^ - M^) 

X [S{p^)n{p, m) + ^{-p°)m{p, p)]}, (8) 

with 


^NJL = — m)2p + jCi 

= -{p{i^ - m)ip + g[{'(pipf' + ( 1 ) 

The column vector ^ = {u,d) represents the quark fields, 
where we take the number of flavors Nf = 2, and the 
number of colors Nc = 3 throughout this work, g is the 
effective coupling constant of the four fermion interac¬ 
tion, and T are the Pauli matrices in the SU(2) flavor 
space. The Fierz transformation of jCj is 

^IF + 2{lpij5T1pf - 2{lpT1pf 

—2('0i75r'0)^ — 4('07^'0)^ 

-4(Vii7^75'0)^ + (2) 

We adopt only and '07o'0 as mean fields on account of 
symmetry properties of the vacuum at finite density, and 
ignore other terms. Then the interaction term is written 

Cimf =2G{'tpip)'tpip - ^77 o ^/’)^7 o ^/’ 

-Giijtpf + (3) 

where G = (4Ac + 1)^/4Ac is the renormalized coupling 
constant, and we define the thermal expectation value of 
an operator O by 


Tr O 

'Tj. 

By using Fierz transformation and mean field approx¬ 
imation we get the Hamiltonian density: 

n = -i'tp'y-\/'tp + Mip'ip+p^a2Af + Ga'l--^a2, (5) 

where A = is the quark number density. 

M is the constituent mass (the effective quark mass): 


M = m — 2Gcri(M, pr). (6) 


a I and <72 are defined as follows: ai = ('0'^), <72 = 
{'ipjo'ip). We also introduce the renormalized chemical 
potential: 


Pr = P - Pj-(72{M,Pr)- (7) 

To calculate M, pr self-consistently by Eqs. (6) and 
(7) we use the formalism of the thermal Green function in 


n{p, p) 


1 

1 qI3{E-h) ’ 


m{p,p) 


1 

1 -h 


(9) 


(3 = 1/T, and E = Using the Green function 

we get (7i and <72 as follows: 


NrNf r V 

cri= -M -^ / —[1 -n(p,/ic) - m{p, gr)]dp,(10) 

Jo ^ 

fA 

/ p‘^[n{p,pr)-m{p,pr)]dp. (11) 

Jo 


NcNf 


In this paper, we employ the parameter set used by 
Hastuda and Kunihiro [18]: m = 5.5 MeV, g = 5.074 x 
10“^ MeV“^, A = 631 MeV, which are determined by 
fitting the pion mass = 138 MeV, the pion decay 
constant A = 93.1 MeV. For the quark condensates we 

obtain = —331 MeV. 

Now one can get the effective quark mass and the quark 
number density for each temperature and chemical po¬ 
tential by solving self-consistent Eqs. (6), (7), (10), and 
(11). In the multi-solution region of Nambu solution (chi¬ 
ral symmetry broken) and Wigner solution (chiral sym¬ 
metry partially restored), the solution which minimizes 
the thermodynamical potential density uj is stable. Erom 
Eq. (5) the thermodynamical potential density uj is given 
by 


U) = A iog[Tr 

9 2 NcNf 


f 


p^Edp 


A.A. 


fr 


TT^ 

Tog[l 


pA 

' / p^{log[l + 

Jo 






( 12 ) 


Having taken this into consideration, we obtain the ef¬ 
fective quark mass in the T — g plane in Eig. 1 in the 
massive and massless current quark case, respectively. 


III. THE CHIRAL PHASE TRANSITION 

The phase diagram of QGD within the two-flavor NJL 
model is shown in Eig. 2. It should be noticed that 
there are multi-solution around the first-order transition 
line [19-24]. The multi-solution of two solutions indi¬ 
cates the competition of the two corresponding phases. 
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FIG. 1. (color online). The effective quark mass M in the 
T — jjL plane. The upper panel corresponds to the case of 
finite quark mass and the lower panel to the chiral limit. 



FIG. 2. (color online). Phase diagram in the two-flavor NJL 
model. The green (blue in the chiral limit) points represent 
the first-order phase transition, the pink points the crossover 
transition, and the red points the second-order phase transi¬ 
tion. 


which is a characteristic of the first-order phase tran¬ 
sition. Such a consideration is very important for de¬ 
termining the first-order transition line and locating the 
CEP. The boundaries of multi-solution region have been 
determined in the Refs. [25, 26] by the approach of NJL 
model and the Dyson-Schwinger equations, respectively. 

We display the multi-solution region in Fig. 3, for mas¬ 
sive and massless current quark, respectively. For the 
massive current quark case, the borders of the multi¬ 
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FIG. 3. (color online). Phase diagram in the temperature- 
chemical potential plane for strongly interacting quarks 
around the GEP/TGP. The upper panel corresponds to the 
case of finite quark mass and the lower panel to the chiral 
limit. 


solution region are marked by the blue and red lines in 
Eig. 3. The blue line is where the Wigner solution begins 
to appear, while the red one is where the Nambu solution 
disappears. The multi-solution region is split into two 
parts by the first-order transition line (the green line in 
Eig. 3). On this line, the thermodynamical potential den¬ 
sities of the Nambu phase {ujn) and Wigner phase {ujw) 
are degenerate. The corresponding Nambu phase is sta¬ 
ble and the Wigner phase is metastable in the left part, 
while the Wigner phase is stabler in the right than the 
Nambu one. With T increasing, the multi-solution region 
tapers off on approaching the CEP (T^^^ ^ 32 MeV, 
//CEP _ 347 MeV). 

Eor the massless current quark case, the massless quark 
is the trivial solution for each temperature and chemical 
potential. The first-order phase transition ends at the 
TCP (T^CP ^ 73 MeV, - 300 MeV), In the do¬ 
main where T > and /i < a second-order 

phase transition instead of the crossover transition oc¬ 
curs. 
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IV. SUSCEPTIBILITIES AND CRITICAL 
EXPONENTS 

A. Susceptibilities 

As the linear response of the physical system to some 
external field, susceptibility is often measured to study 
the properties of the related system. Therefore the stud¬ 
ies of various susceptibilities are very important on the 
theoretical side, which are widely used to study the phase 
transitions of strongly interacting matter [27]. So here we 
introduce several kinds of susceptibilities: the chiral sus¬ 
ceptibility Xs, the quark number susceptibility Xqi fhe 
thermal susceptibility xt, the vector-scalar susceptibil¬ 
ity Xvs : the susceptibilities Xm and Xn (for mathematical 
convenience) defined as follows [28]: 


= 

d{ipip) 


dm ’ 

’ 

XT = 

d{'ip'ip) 


dr ’ 

_ 

Xvs - 7 

dfi 

Xm — 



dm 

Xn Qrji 


From the viewpoint of statistical mechanics it is easy to 
find that Xm equals Xvs- Xm = Xvs = where 

Z is the QCD partition function. 

Here we only display the quark number susceptibility 
in the T — ji plane in Fig. 4 as an example. In the case of 
finite current quark mass, all susceptibilities, i.e., Xs, Xq^ 
XTi Xvsi Xm and Xn change discontinuously on the first- 
order phase transition line, and all diverge at the CEP. 
Then the susceptibilities are continuous in the crossover 
region. 

In the chiral limit, there is a second-order phase transi¬ 
tion line instead of the crossover, and the susceptibilities 
always change discontinuously no matter on the first- or 
second-order phase transition line. All those suscepti¬ 
bilities diverge at the same point, where the first-order 
transition line ends. Therefore they all suggest the same 
TCP. The susceptibilities show different features on the 
second-order phase transition line: Xs? Xt and Xm di¬ 
verge, while Xq and Xn remain finite. The different fea¬ 
tures of the various susceptibilities result in the differ¬ 
ences of the corresponding critical exponents calculated 
in Sec. IV B. 


B. Critical exponents 

The susceptibility in the vicinity of the CEP/TCP di¬ 
verges with an index, the so-called critical exponent. The 
values of these exponents are completely dependent on 
the dimension of space and components of the order pa¬ 
rameter instead of the details of the microscopic dynam¬ 
ics. Hence all theories can be categorized into a much 
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FIG. 4. (color online). 3D plot for the quark number suscepti¬ 
bility Xq around the first-order transition line and CEP/TCP 
in the T — fi plane. The upper panel corresponds to the case 
of finite quark mass and the lower panel to the chiral limit. 

smaller number of universality classes [29]. However, the 
critical exponents we have obtained are mean field val¬ 
ues, due to the mean field approximation employed in 
this work. 

The strength of the divergence is dependent on the 
path, along which we approach the critical points [30]. 
We mark the critical exponents by 7 ^^^, 7 ^^^ for the 
massive and massless current quark, respectively. As is 
well known, the critical exponents of quark number sus¬ 
ceptibility are expected 7 ™^ = 7 ™^ = 1 for the path 
asymptotically parallel to the first-order transition line. 
Eor other paths the critical exponents are expected to be 
7 CEP = 2/3 and 7 ™^ = 1/2 [15]. 

Here we calculate four directions of the critical expo¬ 
nents of susceptibilities: the path from lower /i toward fjp 
(represents or with the temperature fixed at 

(represents or T^^^), we mark this path by 

Analogously, other arrows i represent the path 

from higher /i toward , the path from lower and higher 
T toward T^, respectively. Using the linear logarithmic 
fit we obtain 

lux = -yin |T - T^l -h Cl, (14) 

Inx =-qlnlM - 7l+C2, (15) 

for the directions parallel to the // axis and T axis, respec¬ 
tively. Cl, C 2 are constants and 7 is the critical exponent. 
The fitting procedure of the critical exponent for quark 
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FIG. 5. (color online). The logarithm value of the quark 
number susceptibility as a function of log \fi — at the 

fixed temperature 


number susceptibility for the direction is shown in 
Fig. 5. 

For the case of finite current quark mass, the order pa¬ 
rameter only carries one component—sigma (cr), which 
is the only field becomes massless at the CEP. As is 
widely accepted the phase transition falls into the three- 
dimensional (3D) Ising model as the liquid-gas phase 
transition [31-33]. Our results presented in TABLE. I 
show that, the critical exponents of and 7 ™^ agree 
with the mean field prediction 2/3 for all four different 
directions to approach the CEP. As far as we know, there 
is no discussion on the critical exponents of 7 ™^, 7 ^^^ 
and 7 ^^^ in previous literatures, our results imply that 
they are also 2/3. Because of the finite current quark 
mass and quark chemical potential, the critical fluctua¬ 
tion mixes in the quark number density fluctuation and 
chiral condensate fluctuation as well as in the thermal 
fluctuation, and therefore all the relevant susceptibilities 
show the same critical behavior. Eurthermore, the “t” 
and “H” fields of the 3D Ising model are mapped as cer¬ 
tain linear combinations of T and (i on the T — (i plane 
near the CEP. The singular behavior of the susceptibili¬ 
ties is insensitive to the way of approaching the CEP. It 
only depends on the fact whether the approach is tan¬ 
gential or non-tangential to the “t” direction. Therefore 
it is possible that this character holds even beyond the 
mean field approximation (nevertheless, the real scenario 
need to be checked in further study). 

Eor the reason that the critical region around the TCP 
is chopped off in the chirally symmetric phase, we only 
calculated two directions of the critical exponents, i.e., 
^ and listed in TABLE. II in the chiral limit. 

We obtain and 7 ™^ which agree with the pre¬ 

diction 7 j^^ = 1 [34], and 7 ™^ = 1/2 [34, 35] for 
both directions to approach the TCP. We also calcu¬ 
late the critical exponent of the order parameter ('^' 0 ) 
marked by /3, which agrees with the Ginzburg-Landau 
effective theory prediction 1/4 [34]. In the vicinity of 
the TCP ^ | 7 ^tcp _ 7 "|i/ 4 ^ ^g reasonable 


TABLE 1. Critical exponents at the CEP. 


Quantity 

Path 

Numerical result 

MF exponent 



0.69±0.01 



t 

0.65±0.01 

0.69±0.01 

2/3 



0.64±0.01 



—> 

0.67±0.01 



t 

0.67±0.01 

0.67±0.01 

2/3 



0.66±0.01 



—> 

0.64±0.01 


^OEP 

It 


0.69±0.01 

_ 

t 

0.66±0.01 




0.68±0.01 




0.68±0.01 


^CEP 

im 

t 

0.66±0.01 

0.69±0.01 

— 



0.65±0.01 



—> 

0.66±0.01 


^CEP 

In 

t 

0.67±0.01 

0.65±0.01 

— 



0.68±0.01 



to expect the critical exponents of susceptibility xt is 
equal to 3/4, since xt = d{^ 2 p)/dT. Analogously, the 
critical exponents of Xm is also expected to be 3/4, as 
Xm = Xvs = /dfi. Our numerical calculations sug¬ 

gest that the critical exponents of 7 ^^^ and 7 ™^ are 
3/4, while 7 ^^^ = 1/2 . More studies related to the 
critical exponents near CEP and TCP can be found in 
Refs. [13, 36-39]. Whereas the critical exponents of TCP 
are classical in three dimensions (since the upper criti¬ 
cal dimension is 3), this is not the case for the critical 
exponents of the CEP. The critical exponents in the uni¬ 
versality class of Z 2 are substantially different from the 
mean field exponents. Mean field exponents are classical 
for D=4, and do not correspond to any D=3 universality 
class. So it is worth noting that the critical exponents 
obtained in this paper are the mean field ones, due to the 
mean field approximation we used in this work. 

Moveover, in the chiral limit, 0(4) ^ SU{2)vxSU{2)a 
symmetry is restored on the second-order phase transi¬ 
tion line [40]. That is, the order parameter carry four 
components—sigma and three pions (cr, tt), which all be¬ 
come massless at those critical points. The critical expo¬ 
nents of Heisenberg 0(4) model has been expected [41]. 
We choose the point (T^ 193 MeV, = 0 MeV), 

which are on the second-order line and attracts much in¬ 
terests in lattice QCD, and calculate the mean field crit¬ 
ical exponents at it. We display the critical exponents at 
this point in TABLE. III. The critical exponent of order 
parameter [3 agrees with the mean field theory prediction 
1/2 [42], so Xt expected to be 1/2, for the same 
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nite temperature and chemical potential in the cases of 
finite current quark mass as well as the chiral limit. We 
have compared the multi-solution region of Nambu and 

TABLE III. Critical exponents at the critical point with zero 


Quantity 

Path 

Numerical result 

ME exponent 


t 

0.50±0.01 

1/2 

7? 

t 

l.OliO.Ol 

1 

It 

t 

0.51±0.01 

— 


FIG. 6. (color online). Critical exponents as a function of /x 
on the second-order line. 


TABLE 11. Critical exponents at the TCP. 


Quantity 

Path 

Numerical result 

ME exponent 

^TCP 


0.25±0.01 

1/4 

t 

0.25±0.01 

^TCP 

—> 

0.99±0.01 

1 

t 

l.OOiO.Ol 


^TCP 


0.51±0.01 

1/2 

t 

0.51±0.01 

^ TCP 

7t 


0.73±0.01 


t 

0.74±0.01 


TCP 

—> 

0.75±0.01 


7m 

t 

0.76±0.01 


TCP 

—> 

0.48±0.01 


In 

t 

0.49±0.01 



reason mentioned for the critical exponents at the TCP. 
We also calculate the critical exponents of some points on 
the second-order line. The result shown in Fig. 6 demon¬ 
strates that most of the points on the second-order line 
share the same critical exponent. However as the point 
gets close to the TCP along the second-order line, the 
critical exponent of the order parameter decreases con¬ 
tinuously from its mean field value to the one at TCP. In 
other words, there is a crossover of different universality 
classes, although the TCP and the 0(4) line are both of 
second-order and their critical exponents are calculated 
in mean field approximation. 


Wigner ones in the massive and massless quark cases in 
the phase diagram. The analysis of the multi-solution 
region is necessary, which provides us the exact first- 
order phase transition line, and then we can calculate 
the susceptibilities and corresponding critical exponents. 
The physical quantities, i.e., the effective quark mass and 
quark number susceptibility in the T — /i plane are pre¬ 
sented for the massive and massless current quark case, 
respectively. We have used the susceptibilities to locate 
the CEP/TCP and studied the critical behaviors in the 
vicinity of them. This study reveals that all susceptibil¬ 
ities show the same critical behavior at CEP, and they 
are all found to be governed by a common critical ex¬ 
ponent 7=2/3, for four different directions to approach 
the CEP. While in the chiral limit, the various suscep¬ 
tibilities show different features along the second-order 
transition line and so are the corresponding critical ex¬ 
ponents. In addition, the critical exponent for order pa¬ 
rameter equals 1/4 at TCP, while is 1/2 on 

the second-order line. But the critical behavior in the 
vicinity of the points on the second phase transition line, 
belongs to the 0(4) universality class, with the mean 
field critical exponent for the order parameter /3^=l/2. 
The different critical exponents of the TCP and the 0(4) 
line, indicates a change of different universality classes. 
Last but not least, the mean field critical exponents is 
presented in this study, while the discussions beyond the 
mean field are also interesting, which no doubt deserves 
further study. 
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